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Abstract 

We show that if two closed hyperbolic surfaces (not necessarily 
orientable or even connected) have the same Laplace spectrum, then 
for every length they have the same number of orientation-preserving 
geodesies and the same number of orientation-reversing geodesies. Re- 
stricted to orientable surfaces, this result reduces to Huber's theo- 
rem of 1959. Appropriately generalized, it extends to hyperbolic 2- 
orbifolds (possibly disconnected). We give examples showing that it 
fails for disconnected flat 2-orbifolds. 



1 Introduction 

We say that two hyperbolic surfaces (assumed closed but not necessarily 
orientable or even connected) are almost conjugate if their closed geodesies 
match, in the sense that for every length / they have the same number of 
orientation-preserving geodesies and the same number of orientation-reversing 
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geodesies. Looking beyond dimension 2, we say that two hyperbolic d- 
manifolds (assumed closed, but not necessarily orientable or even connected) 
are almost conjugate if their geodesies match with respect to length and 
'twist'. The twist of a geodesic (also called its 'holonomy') is measured by 
the conjugacy class in 0(d — 1) of the action of parallel translation around 
the geodesic. To say that geodesies have matching length and twist amounts 
to saying that the corresponding deck transformations are conjugate under 
the action of the full isometry group of hyperbolic <i-space. 

This usage of the term 'almost conjugate' accords with Sunada's well- 
known definition [12] that subgroups of a finite group are almost conjugate 
if they meet each conjugacy class equally. Here, two connected hyperbolic 
surfaces are almost conjugate just if their deck groups intersect equally every 
conjugacy class of the full group of isometries of the hyperbolic plane. Strictly 
speaking, it is the deck groups that are almost conjugate, but it is convenient 
to apply the term abusively to the surfaces themselves — the more so since 
we want to consider disconnected surfaces. 

While we haven't specifically required it in the definition, the matching of 
geodesies between almost conjugate hyperbolic manifolds, whether surfaces 
or manifolds of higher dimension, can and should be taken to respect the 
imprimitivity index of the geodesies as well as their length and twist. 

Please note that here and throughout, by 'geodesies' we mean oriented 
closed geodesies. Because our geodesies carry a designated orientation, the 
number of geodesies of length I will always be even, with each unoriented 
geodesic being counted twice, once for each orientation. So when we say, for 
example, that the number of geodesies of length at most I is asymptotically y, 
we're talking about oriented geodesies; the asymptotic number of unoriented 
geodesies would be |y. 

According to the Selberg trace formula, almost conjugate hyperbolic man- 
ifolds are isospectral: They have the same Laplace eigenvalues with the same 
multiplicity. (Cf. Randol's chapter in [3]; Gangolli [9]; Berard-Bergery pQ.) 
Such manifolds can be constructed using the well-known method of Sunada 
[T2] . Sunada's method is very flexible and powerful, and works in contexts 
that go well beyond the kind of locally homogeneous and isotropic case we 
encounter here in discussing hyperbolic manifolds. It readily yields examples 
of non-isometric hyperbolic manifolds of dimension 2 and higher that are al- 
most conjugate and hence isospectral: For an exposition, see Buser [2J. The 
examples Buser describes are all orientable, but with trivial modifications 
the constructions can be made to yield isospectral pairs of non-orientable 
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surfaces. Since we don't know of any handy reference for this, we'll elaborate 
on this in Section 

We've said that according to the Selberg trace formula, almost conjugate 
manifolds are isospectral. All the examples coming from Sunada's method 
are automatically almost conjugate. This raises the obvious question: 

Question 1 If two hyperbolic manifolds are isospectral, must they be almost 
conjugate ? 

In the case of orientable surfaces, where there is no twisting to contend 
with, the answer is yes: This is 'Huber's Theorem' |10j . dating back to 1959. 
Nowadays we recognize this as a direct consequence of the Selberg trace 
formula, which shows how to 'read off' from the spectrum the lengths of the 
geodesies (and, it should go without saying, the associated multiplicities). 

The purpose of this paper is to prove that the answer is still yes for 
surfaces, even without the orientability assumption. 

Theorem 1 // two hyperbolic surfaces (not necessarily orientable or even 
connected) are isospectral, then they are almost conjugate. 

Here, in contrast to the orientable case, we cannot use the Selberg trace 
formula to read off the lengths of geodesies directly from the spectrum. In 
fact, as we will see, there are 'scenarios' for constructing counter-examples 
consistent with the Selberg trace formula. But the Prime Geodesic Theorem 
comes to our rescue, because we can show that any scenario for constructing 
a counter-example requires the frequent participation of a large number of 
geodesies of length exactly I (specifically, at least a constant times j), and 
having this many geodesies of the same length is forbidden by the Prime 
Geodesic Theorem. 

In higher dimensions Question [1] remains open, even in the case of con- 
nected orientable manifolds. The issues at stake in higher dimensions are 
well illustrated in the proof of Theorem [T] — so you might be interested in this 
theorem even if you don't see why anyone would care about non-orientable 
surfaces. 

To see that the possible existence of isospectral hyperbolic manifolds that 
are not almost conjugate is a question that must be taken seriously, we note 
that among flat manifolds, in dimension d > 3 there exist isospectral pairs 
that are not almost conjugate. Our favorite example of this is the 3-manifold 
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pair 'Tetra and Didi' [5]. In the flat case, some care is needed in defining 
almost conjugacy, because while in a hyperbolic manifold geodesies come 
only in isolation, in a flat manifold geodesies come in parallel families of 
varying dimension. So in the flat case, matching geodesies between manifolds 
involves measuring, rather than just counting. But Tetra and Didi will fail 
to be almost conjugate according to any definition. 

Note. For further insight into the possible existence of isospectral spaces 
that are not almost conjugate, it is natural to expand the class of spaces 
we're considering from manifolds to orbifolds. (Cf. Dryden [8]; Dryden and 
Strohmaier [TJ; Dryden, Greenwald, and Gordon [6]) Of course we need to 
extend the definition of 'almost conjugacy' appropriately. We don't propose 
to discuss orbifolds in detail here, but for the benefit of those familiar with 
orbifolds, we have appended some comments in Section [6] below. Briefly, 
what we find is this: Theorem [1] extends to rule out examples among hyper- 
bolic 2-orbifolds. However, there are examples of isospectral flat 2-orbifolds 
(necessarily disconnected) that are not almost conjugate. And we still don't 
know what happens in the hyperbolic case in dimension > 3. 



2 Outline 

Let M be a hyperbolic surface. Given a geodesic 7 of M, let £(7) denote its 
length and 1/(7) its imprimitivity index (the number of times 7 runs around 
a primitive ancestor). Define the weight wt(7) as follows: 



- if 7 is orientation-preserving; 

wt M (7)H " (!) 
-tanh(Z/2) if 7 is orientation-reversing. 

Define the total weight function 

W M (l) = £ wt( 7 ). (2) 

l(7)=i 

From the work of Berard-Bergery pQ and Gangolli [9J we extract the 
following: 

Proposition 1 Let M and N be hyperbolic surfaces, possibly non-orientable 
or disconnected. M and N are isospectral if and only ifWu = Wjv- 
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We outline the proof in section [3J 

In light of Proposition [TJ to prove Theorem [T] above it will suffice to prove 

Theorem 2 If M and N are hyperbolic surfaces and Wm = Wn, then M 
and N are almost conjugate. 

Observe that this is a purely geometrical statement: All reference to the 
Laplace spectrum has been laundered through the total weight function. 

To prove Theorem [2J we will analyze how we might engineer agreement 
between Wm and Wn without having total agreement between the geodesies 
of M and N, and show that this is not possible without having infinitely 
many lengths I for which the number of geodesies of length exactly I is at 
least Cj, for C > 0. This will contradict the following Proposition, which is 
a simple consequence of the so-called 'Prime Geodesic Theorem'. 

Proposition 2 For any compact hyperbolic surface, the number of geodesies 
of length exactly I is o(y). 

Proof. According to the Prime Geodesic Theorem (see [H]), for a con- 
nected hyperbolic surface (whether orientable or not) the number F(l) of 
geodesies of length at most I is asymptotic to j. The number /(/) of geodesies 
of length exactly / is given by the jump of F at I: 

f(l) = lim F(s) - lim F(s). (3) 

s— >/+ s— >i— 

But if F is any positive increasing function asymptotic to G, the jumps of F 
are o{G). So /(/) = o( j). This establishes our claim for connected surfaces. 
The extension to the general case is immediate, because the o( j) estimate 
holds separately on each of the finitely many connected components. ■ 



3 Proof of Proposition d 

Proposition [1] is the 2-dimensional case of a general result that applies to 
hyperbolic manifolds of any dimension (cf. Berard-Bergery [1]) and indeed 
to quotients of any rank-1 symmetric space (cf. Gangolli [9J). All these 
results grow from the original work of Huber [10] , who proved Proposition [1] 
in the case of connected orientable hyperbolic surfaces, using what amounts 
to a special case of the Selberg trace formula. Huber's proof applies without 
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essential change here, and in the more general cases alluded to above. Huber's 
paper is really beautiful and well worth reading. For convenience we outline 
here how the proof goes. 

In the case Huber was considering, namely orientable surfaces, what we're 
calling here the 'total weight function' reduces to what Huber calls the 'length 
spectrum' (Langenspektrum). (Please note that this is not what various 
authors nowadays call the 'length spectrum'.) Huber's Satz 7 says orientable 
surfaces with the same total weight function have the same Laplace spectrum 
and the same volume. (Actually he says they have the same genus, which 
in his case is equivalent.) His Satz 8 says surfaces with the same Laplace 
spectrum have the same weight function and the same volume. 

Huber proves these results with the aid of the following 'Dirichlet series': 

In terms of the total weight function Wm, this is 

i 

D M (s) = J2 W M (l)l ^ 2 cosh- I. 

Huber's entire argument goes through unscathed here, once we have suitably 
modified this Dirichlet series to take proper account of the reduced spectral 
contribution of orientation-reversing geodesies, and revised the definition of 
the total weight function in a corresponding way. The appropriate series was 
worked out by Berard-Bergery [TJ for a hyperbolic manifold of any dimension 
d: 

D m{s) = £ ^Q(*(7), A( 7 )) cosh- 1( 7 ), 

V "(7) 

where A (7) £ 0(d — 1) represents the action on the normal bundle of trans- 
lation around 7, and 



Q(l,A) 



det I 



1 A + A tr ^ 



cosh/ 



This Dirichlet series converges when s — a + it, a > d — 1. Later Gangolli 
[9] generalized this even further, to handle the case of any rank-1 symmetric 
space. 
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Specializing to surfaces, for an orientation-preserving geodesic (A = [1]) 
the Q factor is 



cosh / 



y cosh I — 1 J 

while for an orientation-reversing geodesic (A = [—1]) it is 
The ratio of the latter to the former is 

i 

/ cosh I — 1 \ 2 / 
= tanh ■ 



ycoshZ+i; 2' 

which is just the factor incorporated into our definition of the weight. With 
this definition, the Dirichlet series becomes 

In terms of the total weight function, we once more get 



i 



Dm(s)=Y,W m (1)1(—^— [ ) cosh-/. 



i 



Huber's argument, which goes through here just as in the orientable case, 
is that this Dirichlet series is spectrally determined, and when restricted to 
a vertical line s — a + it, a > 1 it represents an almost periodic function 
(superposition of a discrete set of sinusoids), whose Fourier coefficients tell the 
total weight function Wm- So the spectrum determines the total weight (and 
the volume). The converse goes by observing that D M (s) has a meromorphic 
extension, from the poles of which we can read off the Laplace spectrum of 
M, together with the volume. So the total weight determines the spectrum. 
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4 Proof of Theorem 
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Let «m(0 denote the number of primitive orientation-preserving geodesies 
in M of length exactly I, and /3m (0 the number of primitive orientation- 
reversing geodesies. 

Fix two surfaces M and iV with Wm = Wn, and set 

a(l) = a M {l) - a N (l); (4) 

b(l) = p N (l) - M (l). (5) 

Note that, in the second definition, M and N have traded places. 

Our job is to show that a(l) = b(l) = for all I. The condition Wm = Wn 
tells us that 

E Ub= E i6(^)tanh(i/(2ife))+ E 
fe ez+ k ocici even 

Note how on the right-hand side we have had to distinguish between odd and 
even k, since going around an orientation-reversing geodesic an even number 
of times yields an orientation-preserving geodesic. 

This system of constraints on the integer-valued functions a and b has 
solutions which at first blush look like they might permit the construction of 
a counter-example. To get the simplest solutions, fix an integer q > 2, and 
set 

= log q, 

so that 

tanh(^) = t^l. 
V 2 ' q n + l 

Let 

c n = -EM^/i)^, 
n t 

3\n 



where fi is the usual Mobius function: 



(— l) k if n is a product of k distinct primes; 
otherwise. 



We get a solution to our equations by taking a(l) = b{l) = when I is not a 
multiple of l , and setting 

o(M = q- 1, 
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b{k) = q + l, 

a(2l ) = 1, 

6(2Z ) = 0, 
a(nZ ) = &(n/o) = c n , n = 3, 5, 7, . . ., 
a(n/o) = b(nl ) = , n = 4, 6, 8, — 

The problem with these solutions is that they grow too fast: For n odd, 
a(nl ) is asymptotic to q n /n. Our proof of Theorem [2] proceeds by showing 
that this kind of runaway growth is unavoidable. 

Note. When q = 2, for the sequence ci, C2, . . . we get 

2,1,2,3,6,9,18,30,56,99,...; 

when g = 3 we get 

3, 3, 8, 18, 48, 116, 312, 810, 2184, 5880, .... 

Looking these sequences up in Neil Sloane's online encyclopedia of integer 
sequences reveals that c n tells the number of primitive length-n necklaces with 
beads of q colors, when turning the necklace over is not allowed. Or in other 
words, c n tells the number of equivalence classes of the set {0, . . . , q — l} n 
under the action of the cyclic group Z/nZ. When q is a prime power, we 
have the alternative interpretation of c n as the number of irreducible monic 
polynomials of degree n over the field with q elements. This seems suggestive: 
We can't use these solutions in the context of isospectral hyperbolic surfaces, 
but perhaps we could get some mileage out of them in a different context. . . . 
We leave that question for another day, and get back to the proof of Theorem 

12 

Let 

L = {I : a(l) ^ or b(l) ^ 0} (6) 

and 

L = {I E L : / is not a multiple of any other element of L}. (7) 

According to this definition, L is the set of lengths of geodesies where M 
and N exhibit different behavior, and L consists of those lengths which 
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are minimal with respect to the partial order where I ^ m means m = kl, 
k G N + . Every element of L sits above some minimal element, i.e. 

L C L N+. (8) 

Our job is to show that Lq = 0. 

Lemma 1 \Lq\ < oo 

Proof. Suppose I G L . By assumption, Wm(1) = Wn(1). Because I is 
minimal in L, any contributions by imprimitive geodesies to Wm(1) are ex- 
actly matched by contributions to Wjy(l). This means that the contributions 
of primitive geodesies of length I must match: 

a(l) = tanh(Z/2)6(Z). (9) 

Assume for convenience that a(l) > 0, and hence b(l) > a(l). Rewrite the 
equation above: 

b(l) -a(l) =6(Z)(1 -tanh(Z/2)); (10) 
6(Q - a(Z) 

[) l-tanh(Z/2)' 1 j 

When Z is large, 

According to Proposition [2J the total number of geodesies of length exactly 
Z is o(j). Here we have at least something on the order of ~ geodesies of 
length I. This puts an upper bound on Z, and thus forces |L | < oo. ■ 
Let P D dd denote the set of odd primes. 

Lemma 2 For any I G L , only a finite number of the odd prime multiples 
of I are also multiples of an element of L differing from I. Specifically, 

|ZP odd n(Lo-{Z})N + |<|L |-l. (13) 

Proof. If li G Lq, l\ 7^ Z, then |ZP D dd H ZiN + | < 1. (This is a simple fact 
about divisibility: It has nothing special to do with lengths of geodesies!) ■ 
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Now fix any I £ L , and let p be an odd prime that avoids the finite set 
for which pi e (L — {/})N + . Since / is minimal in L, as above we have 



a(l) = tanh(Z/2)6(Z). (14) 

As above, assume for convenience that a(l) > 0, and hence b(l) > a(l). 

By assumption, Wm(1) = Wn(1). The only geodesies that are 'in play' 
at length pi are those of length I or pi: That was the whole point of the 
restriction we've placed on p. So 

a(pl) + -ail) = tanhQoZ/2) (b(pl) + -&(/)] . (15) 
P \ P J 

Let's rework this: 

a(pl) - tanh(pZ/2)&(pZ) = -(tanh(pl/2)b(l) - a(l)); (16) 

a(pl) - b(pl) + b(pl)(l - tanhW/2)) = -(tanh(yl/2)b(l) - a(l)); (17) 

P 

I(tanh(p//2)6(0 - o(/)) + (6(pl) - a(pi)) 

= l-tanh( P //2) • (18) 

Look at the numerator here. For p large, ~(tanh(p£/2)&(Z) — a(l)) is close 

to — cb(l)), and 6(/) — a{l) is a positive integer. And b{pl) — a{pl) is 

always an integer: Not necessarily a positive integer, just some integer. As 
soon as p is larger than 2(6(/) — a(/)), ~(b(l) — a(Z) ) will be a positive fraction 
smaller than 1/2, and adding an integer to it can only increase its absolute 
value. This means that for p large, the smallest the numerator can be in 
absolute value is something like ^(b(l) — a(/)), which is at least ~. 

Meanwhile, the denominator is 1 — tanh(p//2) w 2e~ pl . So b(pl) is bigger 
than something like f^-. This contradicts Proposition [2] — unless Lq is empty! 
So L = 0, and M and N are almost conjugate. ■ 



5 Isospectral nonorientable surfaces 

It is well known that there are many examples of isospectral closed hyper- 
bolic surfaces. The first examples go back to Vigneras [13], who constructed 
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arithmetic examples from quaternion algebras. More recent constructions 
have used Sunada's method. Sunada's method is very flexible, and can pro- 
duce nonorientable examples as easily as orientable examples. But as we 
don't know of a reference for this, we briefly outline the procedure here. For 
necessary backgound, see Buser [2]. 

If you take any pair of Sunada isospectral closed hyperbolic surfaces with- 
out boundary, then they have a common quotient. Now just add what Con- 
way calls a 'cross-handle' to this quotient, i.e., take the connected sum with 
a Klein bottle. Or more generally, take the connected sum with any closed 
non-orientable surface. Put the hyperbolic metric on this new quotient, and 
lift everything (the cross-handles and the metric) back up to the covers. The 
resulting surfaces are isospectral and nonorientable. 

Another way of producing isospectral nonorientable pairs is to change 
some of the gluings in known orientable examples where isospectrality is 
proven using transplantation. To take a specific example, consider the sur- 
faces described by Buser [2], Chapter 11, page 304. If you reinterpret Buser 's 
gluing diagrams (Figures 11.5.1 and 11.5.2) so that the identifications on the 
P geodesies are by translation, you get a non-orientable isospectral pair. The 
f3 identifications now add four cross-handles, rather than four handles. The 
transplantation method proving isospectrality in the orientable case contin- 
ues to work here as well. 

6 Comments on orbifolds 

Here, as promised above, are some brief comments about orbifolds. 

There are three independent examples of isospectral flat (disconnected) 
2-orbifolds that are not almost conjugate, one involving quotients of a square 
torus, and two involving quotients of a hexagonal torus. We describe them 
using Conway's orbifold notation [3j. 

A standard square torus has as 2- and 4-fold quotients a 2222 orbifold and 
a 244 orbifold. If we call the torus Si and the quotients S 2 and S 4 , spectrally 

Si + 2S 4 = 3S 2 , (19) 

i.e., you can't hear the difference between a torus with two 244s, and a trio 
of 2222s. 

A standard hexagonal torus Hi has as 2-, 3-, and 6-fold quotients a 2222 
orbifold H 2 (this is a regular tetrahedron); a 333 orbifold H 3 ; and a 236 
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(21) 



(20) 



2Ht + 3H ( 



#1 + 3H 3 
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4^21 



(22) 

(23) 
(24) 



These examples arise from a careful analysis of the contributions of ro- 
tations of various orders to the spectrum via the Selberg trace formula. To 
explain just how this works would take us too far afield. However, it is pos- 
sible to verify isospectrality in these examples in a direct and elementary 
way by using Fourier series to represent explicitly the eigenfunctions of the 
component orbifolds, and checking that eigenvalues match up. 

Among hyperbolic 2-orbifolds, no such examples exist, whether connected 
or not. This is a corollary of Theorem [H together with the observation 
that, in contrast to the flat case, in the hyperbolic case elliptic elements of 
differing order make distinguishable contributions to the spectrum. Again, 
to go further into detail would take us too far afield. 



We're grateful to Emily Proctor, Carolyn Gordon, and Emily Dryden for 
stimulating discussions, and to Humboldt University (specifically the Sonder- 
forschungsbereich 647 project) for hospitality and support during the writing 
of this paper. We're sorry to have to thank Peter Sarnak for pointing out 
that if a function is approaching a limit, its jumps must be getting smaller 
and smaller. 
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